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ABSTRACT 

Multidisciplinary design optimization (MDO) of a 
full vehicle under the constraints of crashworthiness, 
NVH (Noise, Vibration and Harshness), durability, and 
other performance attributes is one of the imperative 
goals for automotive industry. However, it is often 
infeasible due to the lack of computational resources, 
robust simulation capabilities, and efficient 
optimization methodologies. This paper intends to 
move closer towards that goal by using parallel 
computers for the intensive computation and combining 
different approximations for dissimilar analyses in the 
MDO process. The MDO process presented in this^ 
paper is an extension of the previous work reported by 
Sobieski et al In addition to the roof crush, two full 
vehicle crash modes are added: full frontal impact and 
50% frontal offset crash. Instead of using an adaptive 
polynomial response surface method, this paper 
employs a DOE/RSM method for exploring the design 
space and constructing highly nonlinear crash functions. 
Two MDO strategies are used and results are compared. 
This paper demonstrates that with high performance 
computing, a conventionally intractable real world full 
vehicle multidisciplinary optimization problem 
considering all performance attributes with large 
number of design variables become feasible. 

INTRODUCTION 

Continuous demands on efficient design of vehicle 
safety, NVH, durability, and other attribute performance 
have increasingly emphasized on the analysis of the 
vehicle structural designs as well as occupant restraint 
systems. Numerical computation methods have been 
widely employed for this purpose. The most popular 
and flexible computation method for vehicle design is 
the finite element methods (FEM). Over the past ten 
years, tremendous increase in computer speed and rapid 


evolution and development of theoretically sound, 
robust and efficient FEMs for the simulation of 
nonlinear structural dynamics have advanced computer 
aided vehicle design to the point where the results are 
trusted with a high degree of confidence. 

The resulting surge in super computing is 
revolutionizing the way vehicles are designed. The 
application of crashworthiness optimization to vehicle 
design has drawn significant attention and interest in 
automotive industry over the past few years [1-6]. 
Recently, Yang et al. [5] developed a nonlinear 
response surface based safety optimization and 
robustness process, which has been successfully applied 
to the vehicle safety design. They investigated four 
nonlinear response surface methods for different crash 
modes of large-scale systems as well as occupant 
restraint system. Sobieski et aL [7] investigated the 
multidisciplinary design optimization for a car body 
structure under constraints of NVH and roof crush. This 
study extends the previous work to include two more 
crash modes: full frontal impact and 50% frontal offset 
crash. A 512-cpu SGI Origin 2000 computer is used for 
computation as opposed to a 256-cpu one as in [7]. In 
addition, an alternative optimization strategy is 
investigated to compare the results. 

All the crash simulations are performed on NASA 
AMES SGI Origin 2000 machines: LOMAX (512 
processors, 300 MHZ) and STEGER (256 processors, 
250 MHZ). The nonlinear explicit finite element 
commercial software, RADIOS S, is used for crash 
simulations while MSC/NASTRAN is used to perform 
NVH analyses. The conservative Taylor Series 
approximation [8] is used to approximate the NVH 
performance functions. 

VEHICLE MODELS AND DESIGN TARGETS 

Vehicle safety design is one of the major attributes 
in car product development. The vehicle structure must 
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be designed to absorb enough crash energy through d 
structural deformation and attenuate the impact force to 
a tolerable level when crash events occur. In the real 
world, all crash modes need to be considered 
simultaneously for crash analysis and optimization. In 
this paper however, only full front crash impact, 50% 
frontal offset impact, and roof crush are considered, as 
the main scope of this research is to demonstrate the 
state-of-the-art MDO methodologies with high 
performance computing. In addition to the safety 
attributes, the vehicle NVH performance measures are 
also included in this study. 

Full Frontal Crash Model 

The full front car crash finite element model used 
in this study contains about 100,000 elements. It 
crashes into a rigid 90 degree fixed barrier with the 
speed of 35 MPH. The key safety performance measures 
in the full frontal crash include occupant Head Injury 
Criteria (HIC) and Chest G, which are calculated from 
the MADYMO analysis with importing the crash pulse 
from RADIOSS crash analysis. The MADYMO is a 
commercial multi-body occupant simulation product 
from TNO. The process is shown in Figure 1. 

The full frontal crash is commonly used to design 
and validate the vehicle front structures. Federal Motor 
Vehicle Safety Standards 208 (FMVSS) specifies the 
safety regulations and test configuration. The regulation 
states that the HIC and Chest G injury numbers have to 
be within 1000 and 60g, respectively. The explicit 
finite element dynamic software RADIOSS is used to 
perform crash simulations throughout this study. The 
design targets for the full frontal impact in this study 
are to satisfy both FMVSS 208 regulation and 
corporate guidelines with occupant HIC and Chest G 
targeted within 450 and 45g, respectively. Note that the 
numbers may not be realistic, as they are solely used 
for proving this methodology. Another design 
requirement is the New Car Assessment Program 
(NCAP) star-rating criterion, proposed by the National 
Highway Traffic Safety Administration (NHTSA) in 
1994. The NCAP star-rating criterion is derived from 
the total injury probability criteria combining the 
occupant HIC and Chest G numbers. The total occupant 
probability of severe injury is given by; 

where 
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If the occupant P [ota ) is equal to or less than 10%, it 
is graded as 5-star, which is the highest rating in 
NCAP star-rating, with 50% confidence interval 
calculated from student /-distribution, as shown in 
Figure 2. The finite element model is shown in Figure 
3. 

Roof Crush Model 

Vehicle roof crush is a federal mandatory 
requirement intended to enhance passenger protection 
during a rollover event. The test procedure is defined in 
FMVSS 216. The finite element roof crush model for 
this study is converted from a NVH model, as shown 
in Figure 4. 

The explicit finite element dynamic software 
RADIOSS is used for crush simulation. Unnecessary 
parts in the NVH model are deleted and some missing 
parts are added in the roof crush model, e.g., very 
detailed side doors were added and the glasses are 
refined. The total number of elements for roof crush is 
about 120,000. A 72 inches by 30 inches rectangular 
ram is added to perform the 1 roof crush as specified by 
the FMVSS 216. The longitudinal axis of the ram (see 
Figure 4) is at a forward angle (side view) of 5 degrees 
below the horizontal, and is parallel to the vertical 
plane through the vehicle’s longitudinal centerline. The 
lateral axis is at a lateral outboard angle, in the front 
view projection, of 25 degrees below the horizontal. 
The lower surface is tangential to the surface of the 
vehicle and initial contact point is on the longitudinal 
centerline of the lower surface of the ram and 10 inches 
from the forward most point of the centerline. In roof 
crush simulation, the ram normal speed is set to be 7.5 
MPH. 

As described in the FMVSS 216, the force 
generated by vehicle resistance must be greater than 
5,000 lbs (22,240 N) or 1.5 times the vehicle weight, 
w hich eve r is less, through 5 inches of ram 
displacement. In this study, the roof crush resistant 
force is set to be larger than 6,000 lbs (= 27kN). The 
door thickness and material yield stresses are chosen as 
the design variables. 

50% Frontal Offset Crash Model 

In addition to full frontal crash and roof crush, a 
50% offset frontal impact mode is also considered in 
the optimization process. The vehicle finite element 
model is exactly same as the front crash model. The 
only difference is the barrier. In this model, the vehicle 
crashes into a 90 degree fixed rigid wall with 50% 
offset (Figure 5). The impact velocity is 40 mph. The 
RADIOSS is used for the simulation. The key output 
from the frontal offset impact is the toeboard intrusion. 
The design target for toeboard intrusion is set to be less 
than 10 inches. The design variables used for 50% 
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frontal offset crash is same as those used for full front 
crash. 

NVH Model 

In car product development process, different NVH 
models are used for different purposes so that the 
quality of the NVH is high and the cost is at 
minimum. A car body called Body-In-Prime (BIP) is 
used for this study. The BIP is a trimmed body 
without all the closures (door, hood, deck lid) and other 
sub-systems (steering column, fuel tank, and seats) and 
trim items (carpeting, battery, etc.). A trimmed body 
structure may be thought of as a vehicle without the 
suspension and powertrain sub-systems. The BIP can 
also be thought as the "Body-In-White n with glass. The 
BIP plays an important role in determining the 
dynamic characteristics of the vehicle. 

The BIP normal modes, static bending and static 
torsion analyses were conducted using the 
MSC/NASTRAN. The full scale NVH finite element 
model is shown in Figure 6. The total number of shell 
elements is close to 68,000. The total number of nodes 
is about 69,000. The normal modes were calculated 
under the free-free condition. The static bending 
analysis was conduced with front (yz and z) and rear (xz 
and xyz) shock towers constrained while for the static 
torsion rear shock tower supports (xz and xyz) and a 
mid point of the lower radiator support (z) were 
constrained. The bending stiffness calculated using a 
load applied at the front rocker locations was 4,551 
N/mm while the torsion stiffness calculated using a 
torque applied at the front shock tower locations was 
8726 N-m/Deg. The free-free normal mode analysis 
showed that the overall torsion at 26.5Hz and overall 
bending at 38.9 Hz. 

The torsion frequency for the BIP free-free normal 
mode is set to increase by 5% from 26.5 to 27.8 Hz. 
The upper bounds for static torsion and static bending 
displacements are chosen as 3.4 mm and 0.9 mm, 
respectively, i.e., 10% improvement from the initial 
design. 

MPO PROBLEM 

The multidisciplinary design optimization problem 
is to minimize the total vehicle weight subjected to 
design constraints of NVH and 3 safety crash modes: 
roof crush, lull frontal impact, and 50% frontal offset 
crash. The problem is formulated as following: 


Minimize: Vehicle weight 

Subject to: 

NVH.CQiistraiDls: 

21XHz< 

Static torsion ~ D t 
Static bending ~ D b 
where fj = 3 rd frequency 

D t = displacement at local point = 3.4mm 


Db = displacement at local point = 0.9mm 

Ro o f Crush Co ns train t s: 

Crush distance (D) < 5” 

Critical peak load (P cr ) > 27kN (~6000lb) 

Full Frontal Impact Constraints : 

Head injury criteria < 450 
Chest G < 45g 

P total ^ 10 % (i.e. 5-star NCAP rating) 
where P total = Total probability of severe injury 

5 0% Fronta l Q ffstf C r ash C 9 n?fr^ipt§ : 

Toe board intrusion < 10” 

In the multidisciplinary design optimization 
problem, there are 10 global (system) thickness design 
variables including windshield, roof panel, roof rail, 
roof cross members and pillars. The total number of 
design variables for the NVH model is 19, including 10 
for backlite glasses and sheet metal thickness, 9 for the 
stiffness of connection between the backlite glass and 
structures. The thickness design variables contain floor 
panels, jacking/towing on quarter panel, backlite glass, 
shotgun and radiator support. There are 5 subsystem 
thickness design variables for full frontal and 50% 
frontal offset crash models, namely rails and subframe. 

As for the roof crush, 3 thickness and 7 material yield 
stress local design variables are taken into account for 
consideration. 

MPO PROCEDURE 

The multidisciplinary design optimization 
procedure is based on the previous work by Sobieski et 
al. [7], as shown in Figure 7. The difference in this" 1 
paper is that instead of using an adaptive polynonial i 
response surface approximation, a DOE/RSM (Design * 
of Experiments/Response Surface Method) method is 5 
employed to construct the approximation models for 
crash performance functions. Among the various 
methods for DOE and RSM, the optimal Latin ~ 
Hypercube Sampling method is employed to explore ^ 
the design space uniformly to capture the nonlinear 
behavior of crash functions, while the stepwise 
regression method is used to construct the nonlinear 
response surfaces based on the computer experimental 
points. The NVH responses are approximated by the 
conservative Taylor Series Approximation (TSA) as in 

[7]. 

Two optimization strategies are employed to 
perform the inner loop within the MDO process. The 
first takes advantage of the design senitivity analysis 
capability in MSC/NASTTRAN [7] and the second 
takes advantage of both the sensitivity and the 
optimization capabilities. 

In the first strategy, the NVH sensitivities are 
extracted from the MSC/NASTRAN output and 
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approximations are constructed using TSA, In addition 
to the three crash mode responses approximated by the 
quadratic order of stepwise regression, the MDO 
problem is solved by an SQP optimizer. As the NVH 
model is less expensive to run in this case, the NVH 
design variables are updated and the analysis is repeated 
in MCS/NASTRAN for several times (3 in this study) 
while keeping the crash approximation models 
unchanged. After three NVH inner loops, all design 
variables are updated and reflected on the NVH and 
crash models. Simulations are then performed to - 
conform the results from the first MDO cycle. The 
crash approximation models are updated and then 
continue to perform the next MDO cycle if necessary. 
The move limits of the design variables for NVH 
approximation using the TSA are selected to be 20%. 

The second strategy for the inner loop optimization 
is to take advantage of the design optimization 
capabilities in MSC/NASTRAN. Instead of exporting 
the NVH sensitivities and performing inner 
optimization loop manually for three iterations, the % 
optimization process for each MDO cycle is completed 
entirely in the MSC/NASTRAN by imposing the crash 
performance constraints using the DEQATN and 
DRESP2 cards for the explicit crash equations provided 
from the stepwise regression approximations/ The 
advantage of this strategy is that the inner loop can be 
completed without any human intervention. However, 
the approximated crash functions needs to be inputed 
explicitly. 

Optimal Latin Hypercube Sampling (LHS) 

In this research, the optimal Latin Hypercube 
Sampling method [9] is employed to explore the design 
space for constructing the response surfaces for crash 
models. The LHS is chosen due to the absence of a 
prior knowledge of the parametric form of the model. 
The experimental design is directed to minimize the 
bias part of the Mean Square Error (MSE) by 
distributing the sample points uniformly over the entire 
design region. _ The number of runs in LHS is 
determined by the total number of factors including 
control variables and noisy variables. The minimum' 
number of runs is selected to be 3N in this paper, where 
N is the total number of design variables. 

Stepwise Regression (SR) 

Among the various response surface methods 
(RSM), the stepwise regression method is used to 
construct the response surface functions for its 
simplicity and accuracy for structural crash problems 
[5]. The regression analysis techniques have widely 
been used and a detailed description of stepwise 
regression procedures can be found in [10]. In general, 
the stepwise regression model is constructed recursively 
by adding or deleting the independent predictions one 
at a time. When the model is built up, the procedure is 
known as forward selection. The first step is to choose 


one predictor, which provides the best fit. The second 
independent predictor to be added to the regression 
model is the one that provides the best fit in 
conjunction with the first predictor. Given the other 
predictors already in the model, further optimum 
predictors are then added at each step in a recursive 
fashion. Alternatively, backward elimination can be 
used. After certain predictors have been added in the 
model, the predictors are dropped one at a time. The 
predictor that has the least effect on the fit of the model 
is dropped at the stage. The stepwise regression model 
is built by combining the techniques of forward 
selection with backward elimination. 

In this paper, the second order polynomial for the" 
regression is employed to construct the nonlinear 
response surfaces for all crash modes (full frontal, 50% 
frontal offset, and roof crush) except the vehicle weight, * 
where a linear basis function is selected, 

HIGH PERFORMANCE COMPUTING 

The primary computational cost is in performing 
the RADIOSS finite element analysis for the sample set 
of design points corresponding to each of the 3 crash 
modes in section 2. The optimal Latin Hypercube 
sampling method that is used in generating an initial 
set of design points for frontal and offset frontal crash 
and the multilevel orthogonal arrays are for roof crush 
as in reference [7]. The computational details of the 
number of sample points and the elapsed computational 
times are provided in Table 1 . It is important to note i 
that the RADIOSS analysis for the baseline and sample * 
set (3*N) designs can be performed concurrently on a 
multiprocessor machine thereby reducing the elapsed 
time. The two conformation analyses correspond to the 
verification analyses on the optimal designs obtained 

Most RADIOSS crash simulations were performed 
on NASA Ames LOMAX machine (SGI Origin 2000, 
512 processors, 300 MHZ). Each simulation used 4 
processors. Based on numerical experience, running the 
crash simulation with 4 processors using RADIOSS 
code can achieve an approximate speedup of 3, 
compared to a simulation with a single processor. One 
cycle of the MDO process can be completed in 70 hours f 
or approximately 3 days with the 512-cpu LOMAX 
computer running simultaneously. It may require 938 
days to complete if it were executed on a single 
processor. The 938 days estimate is based on the 
following formula: (46*60+46*65+25*70)*3 = 22,500 
hours = 937.5 days (details see Table 1). In other 
words, a total speedup of 321 (22,500 hrs/70 hrs) can 
be achieved. Note that the two conformation runs are 
not counted, as they can only be executed after the 
previous MDO cycle. In Reference [7], a similar MDO 
problem, involving NVH and roof crush, was able to be 
completed in one day while it required 257 days of 
elapsed computing time for a complete solution on a 
single processor of an Origin 2000 server. Reference 
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[6] reported an even better speedup using a more 
advanced computer and a Message Passing Interface 
(MPI) based version of RADIOSS. 

This study shows that with high performance 
computing, conventional intractable vehicle design 
problems now become feasible. 

NUMERICAL RESULTS 

The initial design is started from an infeasible 
region, as shown in Table 2, The constraints of the 
third mode frequency, torsion and bending 
displacements, HIC and toe-board intrusion are all 
violated. The design variables, their lower bounds and 
upper bounds are in Table 3. 

After two MDO cycles, all constraints are satisfied 
and it is shown that the two MDO strategies yield 
comparable results. The total vehicle weight is reduced 
by 14.8kg and 15.6kg* respectively. The objective and 
the maximum constraint histories are shown in Figure 
8. The design histories for the two strategies are 
summarized in Tables 2. The design variables for both 
cases are in Tables 3. The DOE results for frontal crash 
(Pioiai vs. weight) and offset crash (maximum intrusion 
vs. weight) are shown in Figures 9 and 10, 
respectively. It is observed that significant design 
improvements for both strategies are achieved, after two 
complete cycles. In frontal crash, the P, ot ai (total 
probability of severe injury) is improved from 10% to 
8.0% and 7.5%, respectively. While in the offset crash, 
the maximum toe-board intrusion is successfully 
controlled within 10 inches, while reducing vehicle * 
weight. The roof crush performance constraint is 
insignificant as all designs are feasible before and after 
the optimization process. All NVH targets are also met, 

i.e. improving the torsion and bending stiffness by s 
10% and increasing the third mode frequency by 5%. 

CONCLUSIONS 

This research has successfully demonstrated the 
feasibility and benefits of the multidisciplinary design 
optimization methodology with high performance 
computing. The results showed that the MDO/HPC ^ 
methodology could substantially reduce the design % 
cycle time and vehicle weight in the development and - 
certification of new vehicle designs while satisfying the 
functionality requirements. 

Two optimization strategies were employed to 
perform the MDO process and they produced 
comparable results. Optimal designs were achieved with 
significant reduction of vehicle weight in both 
strategies. Apparently, more performance attributes 
(such as durability) and safety crash modes (such as 
side impact, rear impact etc.) with a larger number of 
design variables need to be incorporated into this *■ 
process in the future to solve a real world problem. 
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Figure 3. Frontal Crash FE Model 


Figure 4. Roof Crush FE Model 
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Figure 5. 50% Frontal Offset Crash Model 


Figure 6. NVH FE Model 
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Figure 9. DOE Matrix for Frontal Crash (P lo t a i vs Weight) 
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Figure 10. DOE Matrix for Offset Crash (Maximum Intrusion vs Weight) 
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Table 1 . Computational Requirements for Safety Disciplines 


Crash Modes 

Number of Design 
Variables (N) 

Number of RADIOSS 
simulations 

Total Elapsed Time 
(using 4 CPU/simulation) 
on Origin 2000 machine 

Frontal Crash 

10 (system) + 
5 (local) 

3*N + 1 (baseline) + 
2 (conformation) = 48 

60 * 48 = 2,880 hrs 

Offset Crash 

10 (system) + 
5 (local) 

3*N + 1 (baseline) + 
2 (conformation) = 48 

65 * 48 = 3,120 hrs 

Roof Crush 

10 (system) + 
10 (local) 

L24 + 1 (baseline) + 
2 (conformation) = 27 

70 * 27 = 1,890 hrs 

Total 



7,890 hrs 


Table 2. MDO Design History 


Attribute 

Performance 

Baseline 

Target 

Original* 
Cycle 1 Cycle 2 

Alternative** 
Cycle 1 Cycle 2 


3 rd frequency (Hz) 

26.5 

27.8 <f 3 < 29.3 

29.3 

29.3 

27.8 

27.8 

NVH 

Torsion disp. 1 (mm) 

3.8 

< 3.4 

3.4 

3.4 

3.4 

3.4 


Torsion disp. 2 (mm) 

-3.8 

> -3.4 

-3.4 

-3.4 

-3.4 

-3.4 


Bending disp. (mm) 

-0.97 

> -0.9 

-0.9 

-0.87 

-0.9 

-0.9 


HIC 



356 

411 

411 

357 

Frontal 

Chest G 

42 

<45 

38 

39 

42 

39 


P total (%) 


£ 10 

7.3 

8.0 

9.2 

7.5 

Roof 

Resistance force (kN) 

34.7 

> 27 

30.5 

31.2 

31.4 

31.3 


Intrusion 1 (in.) 



9.7 

9.9 

10 

9.9 


Intrusion 2 (in.) 



9.4 

10 

9.7 

9.8 

Offset 

Intrusion 3 (in.) 

10.9 

< 10 

9.4 

10 

10 

9.8 


Intrusion 4 (in.) 

10.1 

2 10 

8.9 

9.4 

9.5 

9.4 


Intrusion 5 (in.) 

10.5 

< 10 

9.3 

9.8 

9.9 

9.6 

- - 


T 740.5 ~ 

Minimize.;, 

E2£i 

n^a 

17273 

H’’l 724*9^ 


* 3 iterations for NVH loop 

** 5 iterations for NVH loop in NASTRAN 
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Table 3, MDO Design Variables 


No. 

Design Variable 

Initial 

Lower 

Upper 

Original* 

Alternative** 



Design 

Bound 

Bound 

Cycle 1 

Cycle 2 

Cycle 1 

Cycle 2 



Common Design Variables (10) 



f ! -a 

N V' v; = y >7 

1 

Windshield 

3.8 

2.6 

5.0 

2.64' 

2.64/ 

2.64/ 

2.64/ 

2 

Roof panel 

0.7 

0.6 

1.5 

0.64/ 

0.64/ 

0.64/ 

0.64/ 

3 

Roof rail 

1.0 

0.6 

1.5 

0.64^ 

0.64/ 

1.16* 

0.734/ 

4 

Roof cross member (fr) 

1.0 

0.6 

1.5 

0.64/ 

0.64/ 

0.884/ 

0.64/ 

5 

Roof cross member (rr) 

0.9 

0.6 

1.5 

0.64/ 

0.64/ 

0.74/ 

0.64/ 

6 

A-Pillar 

0.8 

0.6 

1.5 

1. 11* 

0.94* 

0.92* 

0.99* 

7 

B-Pillar 1 

1.0 

0.6 

1.5 

0 . 774 ' 

0.64/ 

1.2* 

1.5* 

8 

B-Pillar 2 

0.8 

0.6 

1.5 

0.614' 

0.64/ 

0.96* 

1.2* 

9 

B-Pillar 3 

1.35 

1.0 

2.0 

1.044^ 

1.04/ 

1.56* 

1.8* 

10 

C-Pillar 

0.8 

0.6 

1.5 

0.64' 

0.64/ 

0.94* 

0.64/ 


" •••'.. . : 

i w : ■ .raftmv-fln a ran 




11 

Rear floor panel 

0.76 

0.5 

1.0 

1.0* 

1.0* 

0.99* 

1.0* 

12 

Rear floor cross 

1.4 

0.8 

2.0 

1.87* 

2.0* 

1.62* 

1.4 


member 








13 

Front floor pan 

0.76 

0.5 

1.0 

0.584/ 

0.54/ 

0.54/ 

0.54/ 

14 

Front floor inner 

1.07 

0.5 

1.5 

1.27* 

1.32* 

1.17* 

1.16* 

15 

Jacking/towing 

0.8 

0.5 

1.5 

1.38* 

1.5* 

0.98* 

1.5* 

16 

Quarter panel 

0.8 

0.5 

1.5 

0.9* 

0.9* 

0.85* 

0.84* 

17 

Backlite glass 

3.8 

2.6 

5.0 

2.64' 

2.64/ 

2.84/ 

2.64/ 

18 

Rear tire cover 

0.75 

0.5 

1.0 

1.0* 

1.0* 

0.86* 

0.9* 

19 

Shotgun 

1.22 

0.9 

1.5 

1.14/ 

1.5* 

1.32* 

0.99* 

20 

Radiator support 

0.76 

0.5 

1.0 

0.584' 

0.54/ 

0.76 

0.54/ 

21 

Top edge (x-comp) 

1073 

750 

1395 

1073 

1073 

1074* 

1074* 

22 

Top edge (y-comp) 

367 

256 

478 

368 

368 

388.9 

401.1* 

23 

Top edge (z-comp) 

2734 

1912 

3554 

2734 

2734 

2734 

2734 

24 

Bottom edge (x-comp) 

1424 

1000 

1850 

1424 

1424 

1424 

1424 

25 

Bottom edge (y-comp) 

487 

340 

630 

487 

487 

487 

487 

26 

Bottom edge (z-comp) 

3628 

2540 

5090 

3628 

3628 

3628 

3628 

27 

Side edges (x-comp) 

1521 

1065 

1977 

1521 

1521 

1521 

1522 

28 

Side edges (y-comp) 

520 

365 

675 

520 

520 

522* 

523* 

29 

Side edges (z-comp) 

3874 

2710 

5035 

3874 

3874 

3874 

3874 

- 


Roof Crush Design Variables (10)'’ 

iiSSSllilfli 



30 

Front door 

0.7 


1.0 

0.44/ 

0.44/ 

MSSMM 

0.44/ 

31 

Front door inner 

0.7 


1.0 

0.44/ 

0.44/ 


0.44/ 

32 

Rear door 



1.3 

0.774/ 

0.74/ 

0.74/ 

0.74/ 

33 

Material for A-Pillar 1 

0.207 

0.192 

0.345 

0.207 

0.207 

0.207 

0.207 

34 

Material for A-Pillar 2 

0.207 

0.192 

0.345 

0.207 

0.207 

0.207 

mm 

35 

Material for A-Pillar 3 

0.207 

0.192 

0.345 

0.207 

0.207 

0.207 

mm 

36 

Material for B-Pillar 1 

0.207 

0.192 

0.345 

0.207 

0.207 

0.207 

mM 

37 

Material for B-Pillar 2 

0.207 

0.192 

0.345 

0.207 

0.207 

0.207 

0.207 

38 

Material for fr. door 1 

0.207 

0.192 

0.345 

0.207 

0.207 

0.207 

mm 

39 

Material for fr. door 2 

0.207 

0.192 

0.345 

0.207 

0.207 

0.207 

0.207 


Frontal Crash and 50% Frontal Offset Crash Design Variables (5) ' 


40 

Subframe 

2.0 

1.0 

3.0 

1.04- 

1.04/ 

1.04/ 

1.04/ ! 

41 

Rail 1 

1.9 

1.0 

3.0 

1.04/ 

1.04/ 

1.04/ 

1.04/ 

42 

Rail 2 

1.9 

1.0 

3.0 

1.04/ 

1.04/ 

1.04/ 

1.04/ 

43 

Rail 3 

1.9 

1.0 

3.0 

1.04/ 

1.04/ 

1.04/ 

1.04/ 

44 

Rail 4 

2.4 

1.0 

3.0 

2.044/ 

2.114/ 

2.054/ 

1.994/ 


* 3 iterations for NVH loop 


** 5 iterations for NVH loop in NASTRAN 
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Prediction of HSPPO Evaluation Cost 


References: 

1. Plassman, G. E., Su, P. S., Computer Sciences Corporation, and Spengler, W. E., COMSO, “A Test 
Function Suite with Benchmark Solutions for Unconstrained Optimization”, December, 2000 

2. Plassman, G. E. and Sreekantamurthy, T., “A Contouring Method for Statistical Interpretation of 
Performance Data for a Stochastic DOE Based Parallel Optimizer”, January, 2001 

Attachments: 

Plots of predicted benchmark cost in elapsed time for serial and ideal parallel exeution on an SGI Origin 
2000 computer, where parallel times presume the number of processors, NP = NX. 

Observations 

The evaluation counts and corresponding elapsed (single processor) times for the suite of benchmark 
functions in appendix B of above Reference 1 were used to predict function evaluation only costs expected 
when evaluating HSPPO according to the plan addressed in above Reference 2. This plan anticipates 
computing the number of function evaluations, NS, given below for the various planned solution of each of 
the problems of the test suite. 

NS=(l+NX)2 5 2 k , k=0,9 

For each benchmark function a total of 25*10*6 solutions would be computed, representing 25 replications 
of problems over 10 levels of NS (k=0,9), where each level is solved by 6 levels of HSPPO iteration, NI. 
Since the value of NS doubles for each higher level, the NS cost for all levels is bounded above by twice 
the value of NS for the highest level. The following predictions are based on this bound. That is the 
predicted cost of all 25*10*6 solutions of each benchmark is based on 25*6*2 times the cost of a single 
solution employing the highest level of NS. 


Cost for all 30 benchmarks by NX and over all NX through 64 (total) 


NX 

2 

4 

8 

16 

32 

64 

Total 

Single 

CPU 

Hours 

1.454172 

2.843939 

10.92589 

64.74165 

456.0479 

3496.68 

4032.693 


Cost for all 23 benchmarks, with Order NX Complexity, by NX and over all NX through 64 (total) 


NX 

2 

4 

8 

16 

32 

64 

Total 

Single 

CPU 

Hours 

0.939142 

1 .72799 

3.724177 

9.652021 

29.73078 

99.10535 

144.8795 


The available budget on the NAS Origin 2000 for this fiscal year is 7500 single CPU hours. Given the 
above estimates, we may want to eliminate general testing of benchmarks whose complexity is nonlinear 
(in our case quadratic) in NX. We could consider running all benchmarks through NX=16. We could also 
consider using fewer data points, that is fewer levels of HSPPO iterations, NI, or levels of allowed function 










evaluations, NS. The addressed evaluation plan anticipates the required NS for effective HSPPO solutions 
will increase in approximate linear fashion^with NX. This may not prove true. ^ ^ 

With our present evaluation plan, we run six tests, with varying levels of HSPPO iterations, for each subset 
of HSPPO tests addressing a fixed combination of benchmark problem, design space dimension, NX, and 
level of allowed function evaluations, NS. This is desirable for at least some tests to determine optimal 
algorithm performance characteristics. However, it may well not be necessary for all problems or all 
problem sizes. Once an optimal NI is determined, or its variation over different benchmarks and size of 
NX understood, the cost of further tests could be reduced by a factor of six (6). 

Note that the average execution cost of benchmarks of linear complexity in NX increases in an 
approximate linear fashion only for relatively small NX. For NX > 8 the cost increase is greater than linear 
in NX. 

Test for problems of dimension NX > 64 will need to be very selective to manage resources. The following 
table provides cost predictions for the cheapest (Hyper-Ellipsoid, #2) and the most expensive 
(Michalewicz, M=10, #21) of the benchmarks which are linear with respect to NX in cost. 


Cost in Single CPU Hours for Full Testing of Two Order NX Complexity Benchmarks for NX up to 1024 


NX or Total 


Hyper-Ellipsoid, #2 (cheap) 


Michalewicz, M-10, #21 (costly) 


0.030786 


0.054047 


0.051484 


0.11369 


0.09468 


0.302844 


16 


0.187218 


0.964163 


32 


0.40406 


3.385943 


64 


0.972902 


12.62922 


128 


2.695735 


48.82744 


256 


8.031792 


192.4867 


512 


27.37219 


760.771 


1024 


99.55476 


3034.576 


Total 


139.3956 


4054.111 






















st in CPU Seconds 


Serial CPU Cost over BM by NX for NR(=25) Evaluations 
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Benchmark Problem (BM) 


•N*=4 


N*=8 NX=16 N>t32 - - -A- - - §£ries6 












Ellapsed Cost in Seconds 



A/ 







Cost in CPU Seconds 


Serial CPU Cost over NX for Selected BM for MX=25) Evaluations 


1.&06 


1.E+05 


1.&04 
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1.E+00 
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Cost in CPU Seconds 


Ideal Parallel (M>=NX) Elapsed Cost over NX for Selected BM for 

M\= 25 ) Evaluations 




